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We theoretically analyze the efficiency of a protocol for creating mesoscopic superpositions of ion 
chains, described in [Phys. Rev. A 84, 063821 (2011)], as a function of the temperature of the 
crystal. The protocol makes use of spin-dependent forces, so that a coherent superposition of the 
electronic states of one ion evolves into an entangled state between the chain's internal and external 
degrees of freedom. Ion Coulomb crystals are well isolated from the external environment, and 
should therefore experience a coherent, unitary evolution, which follows the quench and generates 
structural Schrodinger cat-like states. The temperature of the chain, however, introduces a statistical 
uncertainty in the final state. We characterize the quantum state of the crystal by means of the 
visibility of Ramsey interferometry performed on one ion of the chain, and determine its decay as 
a function of the crystal's initial temperature. This analysis allows one to determine the conditions 
on the chain's initial state in order to efficiently perform the protocol. 

PACS numbers: 03.65.Ud, 42.50. Dv 



I. INTRODUCTION 

The quantum to classical transition is an intriguing 
problem of quantum physics [1] and a central issue of 
quantum-based technologies, where efforts are being in- 
vested in developing protocols for implementing quantum 
dynamics of systems of increasing size [2| ■ Increasing the 
size of a system is usually associated with loss of coher- 
ence: Even when the physical object undergoes unitary 
evolution, the particles composing it can often be seen 
as a reservoir for each individual one [3]. As a result, a 
coherent and localized excitation can dephase on a rate 
which increases with the number of components. Such 
dephasing can be examined using the so-called Loschmidt 
echo @HE]i which can be measured by means of the vis- 
ibility of an interferometric measurement performed on 
the system [7J |S]. 

In crystals of trapped ions such an interferometric mea- 
surement can be performed on an internal transition of 
one ion of the crystal. In Ref. [5] a protocol for imple- 
menting Ramsey interferometry on the ion of the crystal 
has been proposed. It was shown that the visibility of the 
interferometric signal, here corresponding to the occupa- 
tion of one electronic state of the ion, gives information 
on the quantum state of the crystal when analysed as 
a function of the time t elapsed between the two Ram- 
sey pulses. This protocol is at the basis of the proposal 
of Ref. [9] to create a superposition of two different crys- 
talline structures across the linear-zigzag structural tran- 
sition p~0TfT2] : The mesoscopic superposition of a crystal 
in the linear and in the zigzag structure can be accessed 
by driving the electronic transition of one ion of the chain 
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FIG. 1: (Color online) Ramsey interferometry with an ion 
chain whose vibrations are at temperature T. A quench 
across the linear-zigzag instability is performed by exciting 
the central ion with a laser pulse in presence of spin-dependent 
forces. In (a) the collective motion is initially in a thermal 
state of a zigzag structure, and the central ion in the internal 
state \g). A tt/2 laser pulse prepares it in the superposition 
(Is) + I e ))/v2- (b) The ion in state |e) experiences a tighter 
state-dependent potential. The corresponding conditional dy- 
namics entangle the ions' internal and external degrees of 
freedom. A subsequent laser pulse performs a —n/2 rotation 
on the ion's internal transition. The final occupation of the 
ground state \g) as a function of the time t elapsed between 
the two pulses allows one to extract information on quantum 
coherence and entanglement created by these dynamics. 



in a set-up where the trap frequency depends on the elec- 
tronic state [9l[T3j. In these settings, a first laser pulse 
prepares the ion in a coherent superposition of the elec- 
tronic states, which evolves into an entangled state be- 
tween the chain's internal and external degrees of free- 
dom as sketched in Fig. [T] The visibility of the Ramsey 
signal does indeed decay fast with the time elapsed be- 
tween the pulses, even when the chain is prepared at zero 
temperature and the dynamics are purely unitary. How- 
ever, for quenches close to the structural instability the 
visibility exhibits quasiperiodic revivals which are visible 
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at longer times t. These revivals appear with a frequency 
which is determined by the frequency of the zigzag mode 
and persist when the size of the crystal is increased. The 
analysis of the spectrum of the visibility's temporal be- 
haviour shows features that can be associated with the 
presence of entanglement generated by the first quantum 
quench [!%] . 

The studies in Refs. Q3] assume unitary dynamics 
of a chain initially at T — 0. The assumption of uni- 
tary dynamics is reasonable for state-of-the-art ion-trap 
experiments [15H17I , in which the coherence times are of 
the order or larger than the typical time scales required 
for observing the dynamics predicted in Ref. 9J. Much 
more stringent is, however, the condition that the chain 
should be initially at temperature T = 0. A feasibility 
analysis needs in fact a more quantitative statement on 
the temperatures required so that the protocol can be 
successfully performed. The present work extends the 
analysis of Ref. [14! to the case in which the chain is ini- 
tially at finite temperature T, considering temperatures 
that can be achieved by means of Doppler or ground- 
state cooling [18]. The behaviour of the visibility as a 
funtion of T is analysed for small quenches across the 
linear-zigzag phase transition. This allows us to deter- 
mine the experimental requirements on the temperature 
of the chain in order to perform the protocol. 

The article is organized as follows: In Sec. [Iljthe pro- 
posal of Ref. [5] is summarized. The theoretical model is 
presented in Sec. |III| which includes the detailed evalu- 
ation of the visibility signal as a function of the number 
of ions and of the initial temperature of the crystal. The 
behaviour of the visibility is analysed in Sec. |IV| for a 
chain of three ions and different temperatures, and the 
conclusions are drawn in Sec. |V] Theoretical details for 
the derivation of the results in Sec. |III| are given in the 
appendix. 



II. ION COULOMB CRYSTALS IN A 
THERMAL STATE 

In this section we briefly review the physical model and 
the protocol proposed in Ref. [?J]. These are the starting 
points of the analysis performed in the following sections. 



A. Ion Coulomb Crystals in State-dependent Traps 

The system we consider are atomic ions, which are 
confined in an external anisotropic potential. The ions 
have been laser cooled to sufficiently low temperatures, 
so that they perform small vibrations about the equilib- 
rium positions determined by the competition between 
the external trap and the Coulomb repulsion [THl US] • 

In the following we assume N identical ions with 
mass to and charge q, which are confined by a linear Paul 
trap |11] or a Penning trap [21 . The ions are assumed 
to have been prepared at low temperature T by means 



of laser cooling, such that they are aligned along a string 
and perform small vibrations about their equilibrium po- 
sitions 18J. In this work we assume that after laser cool- 
ing all ions are in the same stable electronic state, which 
is denoted by \g). In absence of external perturbations, 
their motion is governed by the Hamiltonian [TU [20] 

3N , v 

This Hamiltonian describes the harmonic vibrations 
about the equilibrium positions along the string, where 
Eq is the energy of the classical ground state, w| are 
the normal mode frequencies of the crystal forming a lin- 
ear (or a zigzag) chain and b^,bj are the corresponding 
bosonic creation and annihilation operators. 

Structural superpositions can be obtained by prepar- 
ing, for instance, one ion in a coherent superposition 
of state \g) and a second stable state |e), in which the 
ion experiences an additional potential due to a dipole 
trap [51 l2"2T[2~i] . When this potential is sufficiently steep, 
the ion's equilibrium position is displaced with respect 
to the case in which all ions are in state \g). Then, 
the long-range Coulomb repulsion causes a distortion of 
the crystalline structure, such that the ground state and 
the normal mode spectrum are now different. Therefore, 
when one ion is in |e), its dynamics are governed by the 
Hamiltonian 

3N f , 

ffe = £ e + £K ■ (2) 

3=1 \ / 

Here, Eq is the energy of the corresponding classical 
ground state, wj are the normal mode frequencies and 
bp,bj are the corresponding bosonic creation and anni- 
hilation operators. By preparing the ion's internal state 
in a linear superposition of \g) and |e), the dynamics gen- 
erated by the state-dependent Hamiltonian entangle the 
electronic and motional degrees of freedom [5] . A probe 
of quantum coherence and entanglement can be obtained 
by performing Ramsey interferometry with a single ion of 
the crystal [5]. The scheme is sketched in the following. 

B. Ramsey Interferometry with Thermal States 

In Ref. [8] it was proposed to use Ramsey interferome- 
try as a tool to probe the dynamics and thermodynamics 
of ion chains close to the zigzag instability. In Ref. [5] 
it was shown that the visibility of the Ramsey signal in 
chains of three ions presents features that can be asso- 
ciated with the creation of a superposition of motional 
states corresponding to linear and zigzag structures. In 
Ref. [14] it was shown that some of these features are 
found for chains of generic size N. In these articles, it 
was assumed that the chain was prepared in the vibra- 
tional ground state, at T — 0. The scope of this paper 



3 



is to analyze how the coherence of the superposition is 
affected when the initial vibrational state is not pure but, 
say, a thermal state, as it is the experimental situation 
after laser cooling the chain. In this subsection we briefly 
review the interferometric scheme and derive the expres- 
sion for the visibility, which is going to be explicitly eval- 
uated in the rest of this work. 

Let the initial state of the ion chain be described by 
the density matrix qq = g(t = 0). This reads 



9o = \g)(g\ ® Po ■ 



(3) 



where 



Po 



Z 6XP 



3l 



is the density matrix for the external degrees of freedom, 
with k,B Boltzmann constant and Z = r Ti{e~ Hg ^ kBT " 1 } 
the partition function. In Fig.JT^i the initial state is taken 
to be a thermally excited zigzag structure. A laser pulse 
applied for a time At drives resonantly the transition 
\g) — s- |e) of the central ion, which we label by jo. As- 
suming that the pulse area corresponds to a 7t/2 rotation 
of the dipole, while its duration At is sufficiently short, so 
that the chain motion can be neglected during Ar, then 
the density matrix immediately after the pulse takes the 
form 



where U = (\e)Je\ + \g) jo (g\ +R*\e) jo (g\ +4 \g) JO (e\)/V2 
is the evolution operator describing the dynamics due to 
the laser pulse and the operator i?k(x) = e lk ' x describes 
the mechanical effect on the crystal associated with the 
absorption of a laser photon. 

The crystal evolves then freely for a time t according 
to the Hamiltonian H = \e)(e\H e + \g)(g\H g , such that 
the density matrix at time t reads g(t) = Q2, with 



02 



E 

P,q=e,g 



P2,i 



and 



P2,ee (t) = -U e (t)R kPo Rtut(t), 



P^,gg 



(*) = \u g {t)pMt) 



-IK 



P2,eg(t) = P { ge (t) = —U e (t)R kPo Ut(t) . 



(4) 
(5) 
(6) 



There, U p (t) = exp(—iH p t/K), while (j>is a phase-shift ap- 
plied when the atom is in the excited state, which allows 
one to perform interferometry. At this stage, atomic mo- 
tion and internal degrees of freedom are entangled by the 
state-dependent evolution. A graphical representation is 
shown in Fig. [T]d, illustrating a coherent superposition 
between zigzag and linear chain. Information about this 
structural superposition can be extracted by measuring 



the probability that the central ion is in the ground state 
after a second pulse has been applied which performs a 
— 7t/2-pulse. The probability that the ground state of the 
ion is occupied after the pulse reads 



b,t) =Tr{ Qf \gUg\} = Tr{p f<gg (t)} . 



(7) 



where Qf is the density matrix after the second pulse and 

Pf-.ggi 1 ) = \ {p2,gg{t) +Rl<P2,ee(t)Rk' 

+ e i<f, p 2 , ge (t)Ru + e-^4,p 2 , eff (t)) , (8) 

with k' the wave vector of the second pulse. Using Eq. ^ 
in Eq. ([7]), the probability can be recast in the form 



M) = -(l + Refect)]) 



(9) 



where 0(t) measures the coherence between ground and 
excited state. It reads 



O(t) = Tr{Rl,U e R k p U g ^ , 



(10) 



and determines the visibility V of the Ramsey signal 
through the relation V = In Sec. Ill we determine 



0(t) as a function of the time elapsed between the two 
pulses and of the initial temperature of the chain, and in 
Sec. IIVI we discuss its behavior as a function of the tem- 
perature T for a chain of three ions in a two-dimensional 
geometry. 



III. EVALUATION OF THE VISIBILITY OF 
THE RAMSEY FRINGES 

In this section we carry out the theoretical evaluation 
of the visibility of the Ramsey signal as a function of 
the temperature T and of the time t elapsed between 
the two Ramsey pulses. The calculation here reported 
extends the one presented in |14j . which was performed 
assuming that the ion crystal is initially prepared in the 
vibrational ground state. We also include the possibility 
of a mechanical effect associated with the absorption and 
emission of a photon of the the pulses. The final result is 
reported in Eq. ( 50 1 . It is valid for a three-dimensional 



geometry, for any number of ions and for any initial tem- 
perature T, as long as the assumption of that the ions 
perform harmonic vibrations about their equilibrium po- 
sitions is valid. 



Some useful relations 



In order to evaluate Eq. (10), we make use of the uni- 



tary transformations which relate the states and normal- 
mode operators of the two structures. These have been 
derived in Ref. 14 , and are reported in this section. 
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We write the ions' positions fj (j 



r' 

the corresponding structures, rf,r|, namely, 



. , 3N) as small 



excursions , away from the equilibrium positions of 



r, = r a . 



■9? = r 7 



'3 - >j I <ij - - j >lh 

where the superscripts <? and e indicate whether the 
central ion is in the ground state \g) or in the excited 
state |e). We denote by d 9 - = r| — r| the equilibrium 
displacements of ion j between the structures. The nor- 
mal modes are obtained by diagonalizing Hamiltonian H s 
(s = g, e), after it has been expanded around the equilib- 
rium positions up to second order in the displacements . 
Thus, they depend on the internal state of the central ion 
and read 



Q? = E M *. 



kllk 



(12) 



where matrix M s is the orthogonal transformation diag- 
onalizing the harmonic part of the potential; details are 
reported in Ref. [2]. Equations (111 and (12) give the 
mapping connecting the normal modes in the two struc- 
tures: 



k 

P 3 = E T J kP k ' 



D 9 

3 



(13) 
(14) 



The vibrational ground state of each structure is denoted 
by |0) s , with s = g,e. They are mapped into one another 
by the transformation [T3J [35] 



\0) g = ZV e (ft,...,f3l N )e A \0) c 
Here, A reads 



(19) 



A = jE^ 6 f & 



(20) 



with Ajk a real and symmetric matrix with elements 



.4 



3 k 



£(0: 



iVlk ■ 



The scalar 



Z = dct 



; 1 - A 



(21) 



(22) 



is warranting the correct normalization, while the term 
T> e (j3f, . . . ,/Jfjv-) is a displacement operator of the 3N 
normal modes, when the structure is the one correspond- 
ing to the central ion being excited. It is defined as 



(23) 



with P? the momentum canonically conjugated to the 
displacement Qf and 



with 



T 3l = E M fej M fei ' 



(15) 



k k 

Correspondingly, the annihilation and creation operators 

bj and bj\ defined by relations 6| = ^Jmuj s - / (2h) [Q s + 

iP s / (muij)] and its adjoint, are related by the Bogoliubov 
transformations 



3 



u ik bl 



k 

zZ Uk ^ b k 

k 



E 



•0}. 



E^ f 



with the real and dimcnsionlcss coefficients 




and displacements /?? = J muj 9 / 2h D a . , such that 



(16a) 
(16b) 

(17a) 
(17b) 



VPm=eKp{ftbt-K'b.)}. 



(24) 



It is useful to introduce the relation between displace- 
ment operators in the basis of normal modes of each 
structure. For a generic displacement A 9 , here given for 
the structure in which the central ion is in the ground 
state, they are related by the equation 



where 



and 



^[A«]=2Im[E 



A| = E(^i+Af%) 



(25) 



(26) 



(27) 



B. Evaluation of the visibility for any initial state 



We now evaluate Eq. (10), whose modulus is the visi- 
bility for an arbitrary initia state. We use that i?k(x) = 
exp(ik • x) is a displacement operator for each normal 
mode, such that i?k(x) = V e {n). Here, 



s , 




(28) 
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where Kj 



k > M h.j) is thc 



\/ ' (/'•• •• 

projection of the wave vector onto the normal mode j, 
assuming that ion jo is illuminated (joa labels the a = 
x,y,z displacement of the ion). 

It is convenient to use a coherent state basis for per- 
forming thc evaluation of Eq. (10). Therefore, we take 



the trace in Eq. ( 10 ) over the basis of coherent states \a) e 
of the harmonic oscillators corresponding to the normal 
modes when the ion is in state e, where \a) e = ®j \oij) , 

such that V { e 3 \aj) |0) e = | ci- 



ties of the trace, we recast Eq. ( 10 ) in the form 
d 6N a 



0(t) 



T 3N 



3/ . Using the cyclic proper- 
j. ( 10 ) in the form 

(29) 



(a\R kPQ Ut RlU e \a) 



The initial density matrix can be expressed in the form 



Pa 



TV 



3JV 



P (X 9 ) \X 9 )(X 9 \ , 



(30) 



(g>,- I A?) is the basis of coherent states 

I I' g 



where \X 9 ) g 

of the harmonic oscillators, corresponding to the nor- 
mal modes when the ion is in state g, and Pq(X 9 ) 
is the Glauber-Sudarshan-P distribution containing the 
information over the initial state 26 , with X 9 = 
(A?,. 



. , X 9 3N ). Using Eq. (30) in Eq. (29), we find 



0(t) 



d eN a fd m X 9 



TV 



3N 



n 



(31) 



e {a\R k \X a ) gg {X a \U$Rl,U e 



This expression contains two matrix elements. We write 
the first one as 

e(a\ R k \\«) g = Ze^ e {a\ V e (K)V e (X e )V e (f3 e )e A |0) e , 

(32) 

where we used \X 9 ) g = V g (X 9 ) |0) . The second matrix 
element in the right-hand side of Eq. (31 ) can be rewrit- 
ten as 



g {X g \Ul Rl,U e 

=Ze -i<p[A«(t)] (0 | 



g (X 9 (t)\Rl \a(t)) e 



Vl{p e )Vl{X e {t))Vl{K')\a{t)) e , 



where 



EN 



:e A < t u kj + X 9 k *e +iu * t 



(33) 
(34) 



Using these results, Eq. (29) can be cast in the form 
d 6N a 



0(t) = 



where 



7T 



:iiV 



d 6N X 9 

—^Z 2 ^P {X 9 ) 



(a\e A ^ 



|e At (*'>Ki)> e 



ip = cp[X 9 ]-<p[X 9 (t)}+<p e -<pe> 



(35) 
(36) 



with 



Im + /3?)Af + K ^t\ > ( 37a ) 

3 3 

Im[^(4 + Pl)Xf (t) + £ , (37b) 

3 3 

and where k' is the displacement due to the emission 



of a photon with wave vector k'. In Eq. (35) we have 
introduced the quantities 



9j = Kj + /3| + Xj , 

e'j = k'j + p° + x^t) 



(38a) 
(38b) 



as well as the operators 

A(8) = \j2 A ^bf - en(W - 0k*) ■ (39) 
jk 

Exchanging the order of thc integrations and evaluating 
the operators, Eq. (35) becomes 



d 6N X 9 

°W = / -^m- Z 2 e^P (X 9 ) 3 a (X 9 ) , (40) 



with 

3a(A 9 ) 



TV 



■3N 



(a e 



A(0) 



3 At(e') 



•■ " ' \a(t)) e 



d m 



a 



TV 



(a\9) ee (6'\a(t)) ee (a\e Am \a)e 



-3N e \ I fee 



x e (a(t)\e AHe ' } Ht)) e 



(41) 



The explicit evaluation of the integral in the variables a 
is reported in Appendix [XJ and leads to the expression 



0{t) 
where 



d 6 ^„,.^^v^ pG . (e0eGWeJs ^- ls 



TV 



■3JY 



Po(X 9 



VdetU 



G(7) = E^«-E h | d 



(42) 
(43) 



jk 



with 7j = 8j, &y Here, fl is a complex symmetric 6AT-by- 
6iV matrix, which reads 



with 



1 - A+ -iA" 
-iA~ 1 + A J 



i("i+"fc)t 



Moreover, s is a 6-/V-dimensional vector given by 



s = 



-iS 1 - 



(44) 



(45) 



(46) 



6 



with 



Sf{6,6'] = S J [e}±S*[9 / }e~ i ^ t 



k 



7j 



(47a) 
(47b) 



Equation ( 42 ) gives the visibility as a function of an arbi- 



trary initial state, for an arbitrary number of ions N and 
accounting for the mechanical effect associated with the 
absorption and emission of a photon of the laser pulse. 



C. Visibility for an initial thermal state 

We now evaluate the visibility when the chain is ini- 
tially in a thermal state, as in Eq. (§); we need to inte- 
grate in Eq. ( [42] ) over the variables X 9 taking the distri- 
bution P (X 9 ) = Uj p o( X j), su ch that [26] 



Wf) 



Tt(nf) 



exp 



with 



(nf) = (bfb°) 



\3L 



1 



(48) 



(49) 



the mean vibrational number of mode b 9 . The integral in 
the variable X 9 is a Gaussian integral and the resulting 
visibility reads: 



0(t) 



Z 2 e^e c exp{i£ T A- 1 £} 
(ni) • • • (n 3N ) Vdct n dot X 



(50) 



This expression is valid for any initial temperature T and 
any number of ions, as long as the harmonic approxima- 
tion at the basis of our model is valid. In Eq. (50) we 



have introduced a series of quantities in order to provide 
a compact form. These quantities are given here in order 
to make the presentation self-consistent. 

The prefactors contain two exponentials, whose expo- 
nents take the form 



$={v{ K ]-ip[K'])/2, 



and 



C = G(C)+G*(C')+E E 

j,k=l a,P=± 



s?[k, K'}[n- 



where 







(51) 



The vector C is conveniently decomposed into three 
parts, 



C=1 + J + JC. 



(52) 



The first term on the right-hand side is given by 

'2^ (lj(C)+Ij(C)e-' lu9 A 
' U ? (C)+/l(C')e +i ^V 



T 



where 



jk 



2 

j 

1 \ - 



(53) 

(54a) 
(54b) 



The second term can be written as 

1 - e" 
v^i ..... 

with 



/3f (e +II V - 1) + ±(J"(«0 - J^(«')e +iw S t ) > 



■^w = y^.( K j( u kj + v kj) ± ^(«fcj - 



(55) 
(56) 



The third term reads 




J 7 k k l 
Jjfc u k I 



-lia/9 v N 
s.j L .„, .„ Jr „ Jjfc I fcZ 

5/[K,«'][n-i]^Yf, 



5? [«,«/] [n- 



where 



Yjl = AjfeUjn, - uij . 

k 



Y 



(57) 



(58) 
(59) 



The matrix X in Eq. ( 50 1 is given by the following 
expression, 




Y° 

1 Ira 




with 



_i -iK-<)t 

_i v° P+ l ("f+ w ™)* 
2 im 

(60) 
(61) 

(62) 

The integration in A 9 is facilitated by changing to real 

j+iyj, thereby introducing 

^££>^_f) (63) 



Y/ m — _ VljYj m , 



and the thermal excitation, 

71m = <^r> 



and imaginary parts of X 9 



, -l 



yxx ~y x y 

yVX yVV 



3 

1 1 

i — i 



and 



C 1 

3 



3 



(64) 



7 



9 


-0.1 


-0.005 





0.02 




v y /(2n) (MHz) 


1.470 


1.545 


1.549 


1.565 




A 


0.005 


0.01 


0.015 


0.02 


0.025 


i/di P /(27t) (kHz 


110 


155 


190 


219 


245 



TABLE I: Conversion table for the dimensionless quantities 
to actual frequencies used for three ions with an axial trap 
frequency of v x — 271 x 1 MHz. The critical frequency is v c = 
y/12/5v x = 2ti x 1.549 MHz. 



IV. RESULTS 

We now analyze the visibility of the Ramsey fringes 
when the central ion is subject to a sequence of two 
Ramsey pulses in presence of a state-dependent poten- 
tial. The results we present are obtained by evaluating 
explicitly the visibility in Eq. ( 50 1 for a given set of pa- 



rameters, assuming that the vibrations along the direc- 
tion perpendicular to the plane of the zigzag are frozen 
out, namely, the motion is effectively confined to the x — y 
plane. We will focus on a chain composed by three ions in 
a linear trap with axial frequency v x and transverse secu- 
lar frequency v y . In the following we will consider values 
of v y close to the critical value v c , separating the linear 
from the zigzag phase [HI [3U] , and use the dimensionless 
parameter 
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(65) 



in order to indicate whether the ions form a linear array 
[g > 0), or a zigzag chain (g < 0). The instability is 
at g = 0. The effect of the quench on the chain, due to 
the internal excitation of the central ion, is represented 
by a shift of the trapping frequency that the central ion 
experiences, denoted by ^di p . The strength of the quench 
is here described by the dimensionless parameter 



A = 



'dip 



(66) 



that is here taken to be positive, A > 0. Hence, when the 
central ion is in the excited state, the trapping potential it 
experiences is steeper. Table [T] reports the experimental 
parameters corresponding to the values of g and A we 
consider in this section. 

The plots we present display the visibility, namely, the 
absolute value of the overlap 0(t) in Eq. (50), as a func- 



tion of the time t elapsed between the two pulses and 
of the temperature. The plots are evaluated for a chain 
composed by three 9 Be + ions, at different values of g and 
A and at different initial temperatures T of the chain. 
In the first part of this section we discard possible me- 
chanical effects of the laser pulse; this situation can be 
realised with suitably tailored excitation schemes, for in- 
stance by taking copropragating laser beams in a Raman 
scheme |15) , or by using radiofrequency fields |27j . In the 
last part we then consider a pulsed excitation in which 



the mechanical effect is relevant and analyse its effect 
over the visibility signal. 

Before we start, some considerations on the choice of 
the parameter A are in order. We first note that the 
model we consider, a crystalline structure where the ions 
perform harmonic vibrations about the equilibrium con- 
ditions, require that anharmonicities are not relevant for 
the dynamics we investigate. This sets in general an up- 
per bound to the choice of the quench's amplitude A. In 
addition, anharmonic corrections are naturally relevant 
very close to the linear-zigzag instability [T^] , so that the 
initial and final state should be sufficiently distant from 
the critical point. Hence, this sets a lower bound to A 
when the quench is performed across the linear-zigzag 
instability, such that the initial state is, say, a zigzag 
structure and the excited state is a linear array. The pa- 
rameters we choose are chosen in accordance with these 
conditions. 



A. Initial thermal excitation 

We assume that the initial state of the crystal is a 
thermal state of the corresponding equilibrium structure 
at a given temperature T. Table [TT] reports the mean 
vibrational number of each normal mode for the values 
of g and T we consider in this section. In the following 
we will see that one normal mode will become important 
in our discussion. For the crystal being in the linear 
structure, this mode is (for the parameters considered) 
the zigzag mode [8| 1 1 2] . Its frequency and its motional 
pattern are displayed in Table [IT] in the top row in the 
upper block. The frequency of the zigzag mode goes to 
zero when approaching the linear-zigzag transition, and 
when the mode crosses the transition it becomes mixed 
with a second normal mode. The motional pattern is 
displayed for two values in the zigzag in the top row of 
the lower two blocks of Table [III For convenience we will 
name in the following also this mode as the zigzag mode 
when the crystal is below the transition. We also will 
also use the term soft mode for this mode. 

The corresponding visibility is displayed in Fig. [2] when 
the strength of the quench is A = 0.025. Panel (a) dis- 
plays the visibility at different temperatures for g = 0.02, 
namely, when the initial and final state of the quench 
correspond to excitations of a linear structure. The visi- 
bility for T — is given by the black line and it oscillates 
between unity and a value above 0.95. The oscillation 
is at the frequency of the zigzag eigenmode, which is ex- 
cited by the quench j!4j . As the temperature is increased 
the visibility decays, it still exhibits a modulation, which 
is markedly at a smaller frequency but at a larger am- 
plitude than in the case at T = 0. The correspond- 
ing maxima are at a time-scale which is independent 
of the temperature and exhibit a double-peak structure, 
which becomes evident at sufficiently large temperatures. 
Panel (b) displays the visibility when the chain is initially 
a zigzag structure and the quench is performed across the 
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TABLE II: Mean vibrational number for each normal mode 
of the different initial structures, determined by the choice 
of g. The corresponding temperatures are given in fjK. In 
the second row, the vibrations of the ions for each mode are 
sketched, and are displayed as sums of the eigenmodes of the 
linear configuration. The modes of the zigzag chain are com- 
posed of two normal modes of the linear chain, one of which 
(denoted by the thicker arrows) has the main contribution. 



linear-zigzag instability. The feature characterising the 
behaviour at T = is the rapid decay of the visibility to 
zero, and then the appearance of revivals the period of 
the soft mode. This feature is independent of the num- 
ber of ions |14j . Increasing the temperature leads to a 
decrease of the amplitude of the revivals, as also visible 
in the inset: The amplitude significantly drops already 
at T — 5 fiK. Panel (c) displays the visibility when the 
quench connects two different zigzag structures. Also in 
this case, at T = the visibility rapidly decays and then 
exhibits periodic revivals. Thermal effects lead to a de- 
cay of the amplitude of the revivals. Here, however, the 
signal is not significantly altered at temperatures as low 
as T = 10 /iK, as one can observe in the inset. 

These features can be better understood by analysing 
the spectrum of the signal. In particular, we choose to 
study the spectrum of the logarithmic visibility [33], de- 
fined by 



Sln(Wra) 



1 

T 



dt ln[V(i)] e" 



(67) 



Figure [3] displays the spectra of the logarithmic visibil- 
ity corresponding to the curves in Fig. |J Let us first 
recall the behaviour of the spectra at T = 0. These 
exhibit well defined peaks at the frequency, or at mul- 
tiples, of the soft mode. As the temperature increases 




FIG. 2: (Color online) Visibility signal as a function of the 
time t elapsed between the two Ramsey pulses. The sig- 
nal is evaluated for temperatures (black line) , 5 fiK 
(brown/dark gray line), 10 fiK (green/medium gray line), 
50 fiK (blue/medium light gray line) and 100 /iK (light 
pink/light gray line). The parameters are A = 0.025 and 
(a) g = 0.02, (b) g = -0.005, (c) g = -0.1. The insets dis- 
play a zoom of (b) the first double-peak and of (c) the third 
peak. 



the corresponding peaks are broadened. Moreover, ad- 
ditional peaks appear that are located at the beat fre- 
quency Wbcat = |wf — | between the zigzag eigenmodes 
of the two structures, namely, the equilibrium structure 
when the ion is in the ground state and the one in which 
the ion is in the excited state. The appearance of a peak 
at this beating frequency is due to the fact that the cor- 
responding mode in the initial configuration is thermally 
excited. The number of peaks increases with the temper- 
ature; they appear at multiples of Wbcat- This behaviour 
shows that the eigenmodes which most relevantly con- 
tribute to the overlap integral, and thus to the visibility, 
are the soft modes of the initial and excited structures, 
while the contribution of the other modes is marginal. 
Note that the soft modes are, for the parameters here 
considered, the ones which are at lowest frequency and 
significantly occupied, as one can see from Table |H} 

We now compare the signals obtained for different 
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FIG. 3: (Color online) Spectra for the logarithmic signals, 
Eq. ( |67[ ), for the visibility curves in Fig. [5] The vertical dash- 
dotted (green) lines show the zigzag eigenfrequency cj| (in (a) 
the line is at 2uf). The dashed (orange) line shows the lo- 
cation of frequency tJboat = \u)f — ujf\. The insets display 
a zoom of the low frequency part of the corresponding spec- 
trum, highlighting the peak structure of the spectrum at mul- 
tiples Of OJbcat- 



strength A of the quench at finite temperature. Fig- 
ure [4] displays the visibility as a function of the elapsed 
time t evaluated for different values of A and when the 
chain is initially at temperature T=100 /iK. The case in 
which the two structures are linear is shown in panel (a). 
Here, the peaks arising from initial thermal occupation 
wander to later times for smaller values of A. This 
can be understood by recalling that these peaks are de- 
termined by the beating between the zigzag modes of 
the initial and of the quenched structure: For weaker 
quenches, the difference becomes smaller and the associ- 




(MS) 



FIG. 4: (Color online) Visibility as a function of the time t 
elapsed between the two Ramsey pulses for T = 100 piK and 
A = 0.005 (light pink), A = 0.010 (blue), A = 0.015 (green), 
A = 0.020 (brown), A = 0.025 (black line). In a scale of 
gray, as A increases the line becomes darker. The quench is 
performed for (a) g = 0.02 and (b) g — —0.1, such that the 
initial and quenched equilibrium strcutures are either both 
linear or zigzag chains. 



ated timescale, which is the period of the beating, corre- 
spondingly longer. Panel (b) shows the behaviour when 
the initial and the quenched structures are both zigzag. 
Here, the visibility rapidly decays to zero, and then ex- 
hibits some revivals whose height also damps down to 
zero. This latter decay is slower for weaker quenches, 
i.e., for smaller values of A. 

The signal at A = 0.005 is singled out in Fig. [5] where 
it is plotted for a larger interval of elapsed times t. In 
Fig. |5a[ where the initial and quenched structures are 
linear, the signal shows a slow modulation and a certain 
regularity. In Fig. |5b[ where both structures are zigzag, 
the main peaks of the signal are less regularly distributed 
and exhibit a fast quasi-periodic modulation (see inset). 
In both situations one observes that at large times the 
visibility can be significantly above zero, showing that 
coherence persists over time scales of the order of mil- 
liseconds. 

The relevant time scales associated with these features 
become evident by studying the spectrum of the loga- 
rithm of the visibility. Figure [6a] displays the spectrum 
corresponding to the signal in Fig. [5a] The additional 
curves, from top to bottom, correspond to decreasing val- 
ues of the temperature. The black solid line reports the 
case T = 0, which is here plotted for comparison. This 
curve displays clear peaks at values of the frequency cor- 
responding to normal modes or to sum or difference of 
normal mode frequencies. The highest peak here corre- 
sponds to the frequency difference |w| — which for 
the linear chain when all ions are in the ground state cor- 
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FIG. 5: (Color online) Visibility as a function of the elapsed 
time t for 100 /iK, A = 0.005, and (a) g = 0.02 and (b) g = 
—0.1. The vertical lines indicate the location offbeat (dashed 
orange) and of the beating frequency between the transverse 
center-of-mass frequencies (dash-dotted green). The inset in 
(b) shows a zoom for the interval centered about ojbeat- 



responds to center-of-mass oscillations in the transverse 
direction. This peak is still present at finite tempera- 
tures but becomes less prominent. On the other hand, 
at finite temperatures one observes the appearance of the 
peak at Wbcat- Moreover, resonances at multiples of Wbeat 
appear and their number increases with the temperature, 
as is evident by inspecting the inset. 

The spectrum in Fig. [6b] corresponds to the case in 
which the quench connects two zigzag configurations. 
Here, one observes that the peaks characterizing the spec- 
trum at T = correspond to the frequencies of the nor- 
mal modes; they are also present at finite temperatures, 
even though they become broader. At finite T a peak ap- 
pears at Wbcat , while the number of harmonics increases 
with T, as visible in the inset. 



B. Discussion 

The evaluated visibility shows that finite temperatures 
lead to the appearance of various features, which emerge 
because of coherence between the two states created by 
the quench. The fact that the initial state is a statistical 
mixture leads to an overall decrease of the entanglement 
created by means of the quench. In particular, already 
at T = 100 fiK several features of the behaviour at zero 
temperature have disappeared. 

One signature of thermal excitation is the appearance 
of peaks at the harmonic of frequency Wbeat- These can be 
suppressed by cooling the zigzag mode, which is majorly 
excited by the quench, to a lower temperature. Figure [7] 
displays the visibility signal as a function of the elapsed 




FIG. 6: (Color online) Spectrum of the logarithm of the 
visibility, Eq. $67\, for A = 0.005 and (a) g = 0.02 (lin- 
ear to linear) and (b) g = —0.1 (zigzag to zigzag). The 
curves correspond to the temperatures T = (black line), 
T — 5 fiK (brown), T = 10 /iK (green), T = 50 /iK (blue) 
and T — 100 /j,K (light pink), corresponding to a scale of grey 
from dark to light. The inset is a zoom in the low-frequency 
part. The vertical dash-dotted (green) line shows the location 
of the zigzag eigenfrequency in panel (b) and of double the 
zigzag frequency in panel (a). The vertical dashed (orange) 
line shows the location of Wheat- 



time when the zigzag mode has been cooled to 10 /iK 
while the other modes are at T = 100 /xK (see the blue 
line) . By comparing this behaviour with the visibility for 
the chain in the thermal state at T = 10 /iK and 100/zK 
we observe in (a) and (b) that for elapsed times of the 
order of tens of microseconds the visibility qualitatively 
reproduces the behaviour found by cooling all modes at 
10 /iK. Thus, for quenches connecting two linear struc- 
tures (case a) or connecting a linear and a zigzag struc- 
ture (case b) the initial excitation of the zigzag mode de- 
termines the visibility behaviour up to times of the order 
of 10 /is. This is also confirmed when comparing with the 
opposite case, in which the whole chain has been cooled 
to 10 /iK except for the zigzag mode, whose vibrational 
excitations follow a thermal distribution corresponding 
to 100 /iK. The visibility in this case behaves similarly 
to the one where all modes are at T = 100 /iK. A dif- 
ferent situation is found when the quench connects two 
zigzag structures and is displayed in panel (c). Here, all 
eigenfrequencies contribute in determining the dynamics 
at low temperatures. 
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FIG. 7: (Color online) Visibility as a function of the elapsed 
time for A = 0.025 , and (a) g = 0.02, (b) g = -0.005, 
(c) g = —0.1. The pink and black lines correspond to the 
crystal initially at T = 10 [iK and 100 /iK, respectively. The 
green (dark gray) line is for the zigzag mode at 100 /iK and 
all other modes cooled at 10 /iK, the blue (medium gray) line 
is for the zigzag mode cooled at 10 /iK and all other modes at 
100 /iK. 



We finally comment on the mechanical effect of light, 
which in certain configurations of laser pulses can also 
contribute to excite normal modes of the crystal. Its ef- 
fect has been extensively studied in Ref. [5] for different 
parameter regimes. Here, we just show how this may 
modify the signals by assuming that the states \g) and 
|e) are two hyperfine states of the ground state multi- 
plet of 9 Bc + which are resonantly driven by two lasers 
of wave vectors k g and k e via a coherent Raman transi- 
tion [15] , The effective wave vector k, which determines 
the mechanical momentum imparted by the light on the 
ion, is the difference between the wave vectors of the two 
beams and can thus range from zero, when the beams are 
copropagating, to twice the wave vector k g . We choose 
that the momentum imparted by the first pulse is equal 
and opposite to the one of the second pulse, k = k', and 
is always along the transverse, y, direction. We consider 
three situations: (i) copropagating beams, i.e., k = 0; 
(ii) orthogonal beams, k g • k e = (with the resulting 
wave vector along the y axis), and (iii) counterpropagat- 
ing beams along y, k = 2k s . In order to single out the 



effect of the photon recoil, we assume that the crystal is 
at temperature T = 0. 

Figure [8a] displays the visibility as a function of the 
elapsed time when the quench connects two linear struc- 
tures. The signal experiences a visible change due to the 
photon recoil, which is about the same magnitude as the 
variation of the signal when there is no photon recoil. 
In particular, one observes the appearance of other fre- 
quencies, which are due to the excitation of other normal 
modes by the pulse. 

The case in which the quench connects two zigzag 
structures is shown in Fig. |8cJ The curves lie on top 
of each other. The effect of the photon recoil is here in- 
significant on the scale of the variation of the signal which 
is originated by the quench. A similar behaviour is en- 
countered when the quench is across the linear-zigzag 
transition, as visible in Fig. |8b[ where the excitation 
due to the quench dominates over the mechanical effect. 
This behaviour is due to the chosen parameters. Smaller 
quenches give rise to signals where the mechanical effects 
become more visible. 



V. CONCLUSIONS 

The dynamical properties of an ion crystal after a 
quench have been theoretically determined, when the 
crystal vibrations are in a thermal state and the quench 
is performed by creating coherent superpositions of mo- 
tional states close to and across the linear-zigzag struc- 
tural transition. These dynamics are partly revealed 
by performing Ramsey interferometry on one ion of the 
chain. The behaviour of the visibility as a function of the 
temperature has been characterised. 

This analysis shows that the dynamics we predict could 
be experimentally observed in small ion chains. The cal- 
culations we performed took parameters which are exper- 
imentally accessible and show that already at the tem- 
peratures typically achieved by means of Doppler cool- 
ing [18], it is possible to measure the features we iden- 
tify and which are related to the existence of mesoscopic 
quantum coherence between the structures. A prerequi- 
site is that unitary evolution is warranted over time scales 
of the order of hundreds of microseconds, which is fre- 
quently fulfilled in state-of-the-art experiments [T5Hl7j . 
Entanglement generated by the quench can be increased 
if the chain has been previously cooled to ultralow tem- 
peratures by means of sideband cooling [HI [18] or EIT 
cooling [55] of the modes of the chain. 

Another important finding is that the dynamics ex- 
cited by the quench is mostly dominated by the zigzag 
mode of the linear chain, which is the mode driving the 
linear-zigzag instability. A check on the Ramsey signal 
suggests that, after Doppler precooling, it would be suf- 
ficient to ground-state cool the zigzag mode, for instance 
by means of sideband cooling, in order to qualitatively 
reproduce the behaviour of the visibility found when the 
chain is at T = 0. 
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FIG. 8: (Color online) Visibility as a function of the time t 
elapsed between the Ramsey pulses for three 9 Be + ions at T = 
OK. The parameters are A = 0.025 and (a) g = 0.02, (b) g = 
—0.005, (c) g = —0.1. A Raman transition of the central ion is 
driven by two laser beams at 313nm which are copropagating 
(black line), orthogonal (green), or counterpropagating (pink 
line). In a scale of gray the three curves range from black to 
light gray. In panels (b,c), the curves practically lie on top of 
each other. 



This study shows that these dynamics can be observed 
in existing experimental setups. Moreover, our formal- 
ism can be directly applied to larger chains and it can 
be extended by taking into account the normal modes 
including the effects of the micromotion ISHj- The 
analysis in this work provides the basis for investigations 
on the onset of thermalization in closed quantum sys- 
tems [311122]. 
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Appendix A: Calculation of the Overlap Integral 

The integral we have to calculate is given by 

d 6N a 
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Evaluating all expressions, we get 
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Merging all terms into a single exponential and sorting 
them by their orders at once, the exponent reads 
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We now express the integration variables by their real 
and imaginary parts, otj = Uj + ivj and a* = Uj — iVj. 
The quadratic term is rewritten as 
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and the linear term as 
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We defined here the complex symmetric matrices Aj k and The result of the integral is given by 
the vectors 5 ■ by 



A% = - {A jk ± A jk ) , Sf = S j [9] ± S* [e'^P (A5) 

Introducing the vector w = (u, v) T where u — 
(ui, . . . , W3jv)j v = ( v ii ■ ■ ■ j V3n), we can write the in- 
tegral as 
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and the integral in the a's, Eq. (41 1, reads 
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For the demonstration of the convergence of the integral 
( A7) we refer the reader to Ref. |14j , where a proof has been 
reported. 
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